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Ph' Abstract 

We calculate direct CP-violating rate asymmetries in charged B — ► PP and B — > VP 
j> \ decays arising from the interference of amplitudes with different strong and CKM phases. 

The perturbative strong phases develop at order a s from absorptive parts of one-loop 
5_i ■ matrix elements of the next-to-leading logarithm corrected effective Hamiltonian. CPT 



constraints are maintained. Based on this model, we find that partial rate asymmetries 
between charge conjugate B^ decays can be as high as 20% for certain channels with 
branching ratios in the 10~ 6 range. Because the cc threshold lies so close to the physical 
momentum scale, the asymmetries depend sensitively on the model assumptions used to 
evaluate the imaginary parts of the matrix elements, in particular, on the internal mo- 
mentum transfer. The charge asymmetries of partial rates would provide unambiguous 
evidence for direct CP violation. 
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1 Introduction 



So far CP violation JT| has been detected only in processes related to K° — K° mixing 
but considerable efforts are being made to investigate it in B decays. While the most 
promising proposal for observing CP violation in the B system involves the mixing 
between neutral B mesons ||, the particular interest in decays of charged B mesons lies 
in their possibilities for establishing the detailed nature of CP violation. Since charged B 
mesons can not mix, a measurement of a CP violating observable in these decays would 
be a clear sign for direct CP violation, which has possibly been found in K decays where 
the measurements of e'/e now seem to favour a non-zero value |4j which is consistent 
with expectations from the standard model and a top quark mass around 150 GeV 

In non-leptonic charged B decays two main categories of direct CP-violating observ- 
ables can be investigated: First, rate asymmetries || 0], 

_ T(B~ -> f) - T(B+ -> f) 
CP T(B~ ->f) + T(B+ /) ' 1 ' 

and second, azimuthal angular correlations || |9| . The latter involve the decay of the B 
meson into two vector particles B — > V1V2 with subsequent decays of V% and V2 || ||. 
By analyzing the azimuthal dependence of the vector meson decay products one can 
then isolate CP odd quantities. The advantage of this method is that the CP violating 
terms occur even when there are no strong phase differences between the interfering 
weak amplitudes. On the other hand the azimuthal correlations can also be present 
when the CP-violating weak phase differences vanish. By measuring their coefficients in 
charge conjugate B^ 1 decays one has the possibility to disentangle the effects of strong 
and weak phases |], . 

The rate asymmetries occur even for final states with spinless particles but require 
both weak and strong phase differences in interfering amplitudes. The weak phase dif- 
ferences arise from the superposition of amplitudes from various penguin diagrams and 
the usual W-exchange (if contributing). It is clear that a significant contribution of pen- 
guin diagrams, and hence of the CKM [|11|] phase differences, is an exceptional case and 
requires either the absence or a strong CKM suppression of the tree contributions (as 
e.g. in charmless b — > s transitions). The strong phase is generated by final state inter- 
actions. At the quark level the strong interaction effects can be modeled perturbatively, 
following Bander, Silverman and Soni ||, by the absorptive part of penguin diagrams. 
There may be additional final state interaction effects at the non-perturbative hadronic 
level (soft final state interactions). These are very difficult to estimate, but since the 
mass of the B is far above the usual resonance region, we expect these phase shifts to 
be small. 

The rate asymmetries for exclusive two-body decays into pseudoscalars are usually 
estimated using either the model of Bauer, Stech and Wirbel |L2| (BSW) based on 



wave functions in the infinite momentum frame, or the perturbative methods developed 
by Brodsky et al. The rate asymmetries acp can be quite large (of the order 



acp ~ 0.1) for some of the final states. However, the corresponding branching fractions 
of these decays are quite small, ranging from 10~ 6 (estimates with the BSW model 
|12|) to 10 -7 (estimates with the Brodsky-Lepage model []13| , [14], pl|). The number 
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of B-mesons required to resolve the CP asymmetry experimentally (proportional to 
(oq P x BR)^ 1 ) is therefore at least of the order of 10 8 . 

Angular correlations and rate asymmetries for B — > VV have been investigated in 
some detail in our previous work 1 10(1 - There we considered all possible decay channels 



of charged B mesons into two vector particle final states which are induced either by 
b — > s or b — > d transitions. We found appreciable rate asymmetries for the transitions 
Br -> if*- + w, 5^ -»• if*" + p°, 5- -»• if*" + if* and B~ -> p~ + w . For these de- 
cays the predicted rates are necessarily small but large enough to be seen in upcoming B 
factories. These encouraging results have lead us to consider the simpler decay channels 
of charged B's into two pseudoscalars (PP) and into one vector plus one pseudoscalar 
particle (VP) within the same framework as in our recent work for the two vector final 
state flIU| . In B — > PP and B — > VP only one helicity amplitude contributes and thus 
no angular correlation coefficients are at our disposal for detecting direct CP violation. 
On the other hand rate asymmetries might be easier to measure. 

In order to systematically take into account the 0(a s ) penguin matrix elements, 
we base our treatment on the next-to-leading logarithmic short distance corrections 



evaluated by Buras et al. fL6|| . We include also some pure penguin modes and give 
estimates of their branching ratios. In this part there is overlap with other work which 
we will mention later when we present our results. Briefly, because we have a more 
complete treatment of the one- loop matrix element, including all the CP constraints of 
|7|, |17| and because we have included all 0(a s ) penguins (including so-called 'hairpin' 
diagrams and pseudoscalar diagrams) our results are more complete than earlier work. 

The remainder of this paper is organized as follows. In Sect. 2 we describe the 
effective weak Hamiltonian and the evaluation of the hadronic matrix elements. The 
final results for the branching ratios and rate differences are discussed in Sect. 3. A 
table of formulae for the various decay modes can be found in the appendix. 



2 The effective Hamiltonian 



In the next two subsections we follow closely our earlier work on B — > VV decays JlO 



To make the paper self-contained we shall repeat some information already given in |TD 



In subsection 2.3 we describe the evaluation of the hadronic matrix elements which are 
relevant for PP and VP final states. 



2.1 Short distance QCD corrections 

For calculations of CP-violating observables it is most convenient to exploit the unitar- 
ity of the CKM matrix and split the effective weak Hamiltonian into two pieces, one 
proportional to v u = V u bV* s (or V u bV* d in the case of b — > d transitions) and the other 
one proportional to v c = V c bV* s (or V c bV* d correspondingly), 

n cS = ±% (v u nS + «cwg) • (2) 
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The two terms (q = u, c) 



w2 = £*(a«)-Q 



(3) 



differ only by the quark content of the local operators, and for our purposes it is sufficient 
to consider only the following four-quark operators: 



0[ q) = s a j»Lqp ■ qpipLb a , Of = s a ^Lq a ■ q^Lbp , 

3 = s a ^Lb a -J2q'f3lv L( l'f3, °i = s a j'"Lb -J2q'^ fl Lq' a , 

q' q' l 4 J 

5 = s a ^Lb a ■ tfplnWp , 6 = s^Lbp ■ £ q' p l^Rq' a ■ 

q' q> 

where L and R are the left- and right-handed projection operators. The operators 
O3, . . . , Oq arise from (QCD) penguin diagrams which contribute at order a s to the 
initial values of the coefficients at \i ~ Mw [0, or through operator mixing during the 



)(</) 



renormalization group summation of short distance QCD corrections [jig] . The usual 
tree-level jy-exchange corresponds to 2 (with c 2 (M w ) = 1 + 0(0;;;)). The renormal- 
ization group evolution from \x ~ M w to /i ~ has been evaluated in next-to-leading 



logarithmic (NLL) precision by Buras et al. |]16 |. These authors also demonstrated how 
the 0(a s ) renormalization scheme dependence can be isolated in terms of a matrix Yij 
by writing 



8i 



An 



■'.1 



(5) 



where the coefficients Cj are scheme independent at this order. The numerical values for 

240 ± 90MeV from the compilation 
4.8 GeV are 



= 350 MeV (which is in accord with A 
of G. Altarelli pf), m t = 150 GeV and /1 



Cl 

S 5 



-0.324 
0.017 , 
0.011 , 



(5)_ 
MS 

m b 

c 2 
c 4 
c 6 
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1.151 , 
-0.038 
-0.047 



(6) 



Contributions from the color magnetic moment operator will always be neglected in 
the following, because already its tree level matrix elements are suppressed by a factor 
a s /47r and it cannot provide interesting absorptive parts in the decays considered here. 



2.2 Quark-level matrix elements 

Working consistently at NLL precision, the matrix elements of 7i e fr are to be treated at 
the one-loop level in order to cancel the scheme dependence from the renormalization 
group evolution. The one-loop matrix elements can be rewritten in terms of the tree- 
level matrix elements of the effective operators: 



<«jW$I&> = E*(aO 



<% + Illy (//,.. .) 



'^'\r)(l) IMtree 



(sq'<f\Of\b) 



(7) 



The functions m.ij are determined by the corresponding renormalized one-loop diagrams 
and depend in general on the scale /i, on the quark masses and momenta, and on the 
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renormalization scheme. The various one-loop diagrams can be grouped into two classes: 
vertex- corrections, where a gluon connects two of the outgoing quark lines (fig. la), and 
penguin diagrams, where a quark-antiquark line closes a loop and emits a gluon, which 
itself decays finally into a quark-antiquark pair (fig. lb). 

When expressing the rhs of (|7|) in terms of the renormalization scheme independent 
coefficients q, the effective coefficients multiplying the matrix elements (sq'q'\Oj 9 ^\b) tree 
become 

(mij-ry) . (8) 



„cff 



C 



The renormalization scheme dependence, which is present in and r^, explicitly 
cancels in the combination — [TB| . 



The effective coefficients c\ 2 receive contributions only from vertex- correction dia- 
grams, which will not be included in the following (see discussion at the end of Sec- 
tion 2.3). For a general SU(N) color group the remaining effective coefficients can be 
brought into the following form 

1 



„eff 



„eff 



,-eff 



„cfT 



c 3 



2N4n 
1 a. 



y c t + Cp) + 



^ + 24^ {Ct + Cp) + "' 



1 a s 



c 5 - 



(c t + c p ) H 



2N4tt 
C6+ 2 4^ (Ct + Cp) + ' 



(9) 



where we have separated the contributions c t and c p from the "tree" operators 0\^% and 
from the penguin operators O3...6, respectively. The ellipses denote further contributions 
from vertex- correction diagrams. 

In addition to the contributions from penguin diagrams with insertions of the tree 
operators Oil 



ct = c 2 - 



10 2 mi / \ 

9 3 /r K m AJ 



1 j 



(10) 



where AFi is defined in [ID], we have evaluated the penguin diagrams for the matrix 
elements of the penguin operators: 



c 3 



280 2 m 

h -in 

9 3 /i 



2 „ mi . „ 

2 + - 1 



m 



AFi 



k 2 



+ (c 4 + c 6 



E 



- + -^-AF 1 - ) 

9 3 [i l K, mr-' 



3 j 



:n) 



Note that the coefficients c| ff depend on k 2 and, as we shall see later, on the qq states 
that are included in the sum over the intermediate states. 



2.3 Hadronic matrix elements in the BSW model 

To take into account long distance QCD effects which build up the hadronic final states, 
we follow Bauer, Stech and Wirbel fl2|]: With the help of the factorization hypothesis 
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20| the three-hadron matrix elements are split into vacuum-meson and meson-meson 
matrix elements of the quark currents entering in Oi, . . . ,0$. In addition, OZI sup- 
pressed form factors and annihilation terms are neglected. In the BSW model, the 
meson-meson matrix elements of the currents are evaluated by overlap integrals of the 
corresponding wave functions and the dependence on the momentum transfer (which 
is equal to the mass of the factorized meson) is modeled by a single-pole ansatz. As 
a first approximation, this calculational scheme provides a reasonable method for esti- 
mating the relative size and phase of the tree and penguin terms that give rise to the 
CP-violating signals. It is known that the BSW matrix elements do not describe fully 



the existing experimental information on decays like B — > J/tp + K, J/ip + K* [21]. But 
this should not matter too much for the estimates of CP effects presented here. 



Compared to the transitions B — > VV treated in our earlier work plj the amplitudes 
for B — > PP and B — > PV have a simpler structure because only one helicity state 
contributes. On the other hand there are additional contributions from the (V + A) 
penguin operator 0$: After Fierz reordering and factorization it contributes in terms 
which involve a matrix element of the quark-density operators between a pseudoscalar 
meson and the vacuum 

(PiM 2 |0 6 |fl) = -2j2((Pi\qbL\0}(M 2 \sq R \B) + {M 2 \qb L \Q,)(P 1 \sq R \B)) . (12) 

Using the Dirac equation, the matrix elements entering here can be rewritten in 
terms of those involving usual (V — A) currents, 

(P 1 M 2 \0 6 \B) = R[P 1 ,M 2 ](P 1 M 2 \O a \B) , (13) 

with 

Here, m q \ (Wfi) and m q 2 are the current masses of the (anti-) quark in the mesons Pi 
and M 2 , respectively, and the the upper (lower) sign is for the PP {PV) final state. 
We use the quark masses m u = = 10 MeV, m s = 200 MeV, m c = 1.5 GeV and 
m b = 4.8 GeV. 

Finally, one arrives at the form 

(PiM 2 \H^\B) = Z?(P 1 \f\Q)(M 2 \j ii \B) 

+ 4 9) (M 2 |^|o)(p 1 |j;is) , (is) 

where and j' are the corresponding (neutral or charged) V—A currents. The factor- 
ization coefficients z[ q ^ and are listed in the appendix. In terms of the form factors 
Fo,i and Ao for the current matrix elements defined by BSW | 12| , this yields 

(P 1 P 2 \H i S\B) = Z^\Ml-Ml)f Pl F {Ml)+Z$\Ml-Ml)f P2 F {Ml) , (16) 
(PiValwglfl) = X^Ml M P , M*) (Z{ q) f P A (M 2 P ) + Z^ frF^M*)) . (17) 

Mb is the mass of the decaying B meson and fp (fy) are the decay constants of the 
pseudoscalar (vector) mesons in the final state. 
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Concerning how 1/N terms are treated in the coefficients (see (|9|) and (pip), it is 
well known |22| that this model has problems accounting for the decays with branching 
ratios which are proportional to the combination c\ + c%/N. This is due to the rather 
small absolute value of this particular combination when using the short-distance QCD 
corrected coefficients^. An analogous effect is also known in nonleptonic D decays |TJ| . 
and several authors advocated a modified procedure to evaluate the factorized ampli- 
tudes [12 , 24j[ : There, only terms which are dominant in the 1/N expansion are taken 



into account. Recently there has been much discussion in the literature concerning these 
issues. Some authors have argued that QCD sum rules validate this procedure [p5[| . We 



also choose this leading 1/N approximation in evaluating the matrix elements of the 
weak Hamiltonian and we use the QCD corrected coefficient functions q given above. 

The strong phase shifts are generated in our model only by the absorptive parts (hard 
final state interactions) of the quark-level matrix elements of the effective Hamiltonian. 
Of course, when factorizing the hadronic matrix elements, all information on the crucial 
value of the momentum transfer k 2 of the gluon in the penguin diagram (fig. lb) is lost. 
While there has been an attempt [H| to model a more realistic momentum distribution 



by taking into account the exchange of a hard gluon, we will use here for simplicity only 
a fixed value of k 2 . From simple two body kinematics [^6| or from the investigations in 



ref. H] one expects k to be typically in the range 



,2 



m b < 2 < m 



k 2 £ . (18) 
4 2 v ; 

The results we shall present are sensitive to k 2 in this range because the cc threshold 
lies between these limits. Arguments have been made that the lower limit is a more 



appropriate choice ]T7|]. In this work we follow [Kj and choose the upper limit for our 
numerical presentation in the tables. However, we have studied the asymmetries as a 
function of k 2 and will show the results later in Fig. 2-4. 

While the next-to-leading logarithmic precision of the effective Hamiltonian allows 
one to consistently calculate all amplitudes at order a s and to include all one-loop 
matrix elements, some care is necessary when evaluating CP-violating asymmetries of 
the decay rates. In particular, one should make sure that the rate asymmetries for 
sufficiently inclusive channels remain consistent with CPT constraints in certain mass 



limits 0. This issue has been discussed in some detail in |10[ and we shall follow a 
similar procedure of neglecting absorptive parts from flavor diagonal rescattering. 

We refined here the approximation of ref. |10| by dropping only those fractions of 
the imaginary parts of AFi(k 2 /m 2 ) in ( |i~0|) and ( |TT| ) which correspond to the projection 
of the intermediate quark state containing a gg-pair onto the final state. Moreover, no 
absorptive parts of vertex correction diagrams are included in our calculations, because 
they are always flavour diagonal. Of course, this approach may overestimate the CPT 
cancellation of the absorptive phases. We have explicitly checked for all decay channels 
that taking into account the full imaginary parts of the penguin-like matrix elements 
(|10|) and (0) — including flavour-diagonal contributions — does not significantly change 

3 An interesting analysis of the value of c\ + C2/N has recently been presented in ref. |^3| pointing 
out the renormalization dependence of the QCD short distance corrections 
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our results. For the most sensitive cases B — > n^r]' and 5 — > p~r( the inclusion of all 
intermediate flavours would shift the asymmetries by at most —10%, which is in these 
cases still less than the uncertainties from the value of k 2 . 

As a result of our procedure, the factorization coefficients z'fl depend not only on 
whether the basic tree level operator is of the uu or cc type, and on the effective value 
of k 2 , but also on the quark content (uu, dd, ss, cc) contributing in the one-loop penguin 
matrix elements. As in we do not explicitly drop higher order terms which arise, for 
instance, through interferences among (real and imaginary parts of) the order a s matrix 
elements. However, such terms can not introduce the above mentioned inconsistencies 
with CPT because the flavour-diagonal absorptive parts are discarded. 



3 Results and Discussion 

For a numerical analysis of the decay parameters and their CP-violating effects within 
our model, we need to specify the CKM matrix elements and the current form factors. 
It is well known [27] that fits for the parameters Q 

p = COs5 13 S 13 /(s 12 S 23) 

V = sin5 13 s 13 /(s 12 s 23) (19) 

of the CKM matrix depend critically on the value of the 5-meson decay constant f B . 
The solution for lower f B values leads to a negative p while higher f B values render p 
positive. 

We have calculated our results for the two solutions, with the values 

p = 0.32, r/ = 0.31, (f B = 250 Me V, m t = 135 ±27 Gel/) 
p = -0.41 , r] = 0.18, (f B = 125 Me V, m t = 172 ± 15 GeV) 



from the analysis by Schmidtler and Schubert ||27||. A more recent analysis by Ali and 



London based on the latest information on V u b yields similar results. 

The results for the modes which get contributions from tree and penguin diagrams 
are collected in Tab. 1, 2, 3 and 4. They are obtained with k 2 = m 2 /2 as in our earlier 
work [|IIJ and for the two CKM parametrizations which we call the p > and p < 
solution, in the following discussion. The results for B — > PP are given in Tab. 1 
(b — > s transitions) and Tab. 2 (b — *■ d transitions). B — > PV results are given in 
Tab. 3 (b —>■ s transitions) and Tab. 4 (b — > d transitions). The pattern of branching 
ratios of the various channels is rather similar for B — > PP, B — > PV and B — > VV 



Compared to B — > VV decays considered in [T(| the branching ratios are similar but the 



asymmetries are somewhat smaller. Nevertheless there are some interesting cases with 
large enough branching ratios and asymmetries to be measurable in currently planned 
machines. Some differences between the pattern for B — > PP and VP and that for 
B — > VV are due to the contribution of the operator 6 (see (|T3"D). 

If we look at specific decay channels the most interesting cases in Tab. 1 are the 
decays B~ — > K~ + v( and B~ — > K~ + 7r°. The branching ratios are of the order 



These coincide with the parameters p and r\ of the Wolfenstein representation for small angles 
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of 10~ 5 and the asymmetries in the several percent range. For p > they are larger 
than for p < whereas the branching ratios show the opposite pattern. The channel 
B~ — > 7r~ + K° is a pure penguin transition with a moderately large branching ratio 
but a small asymmetry. In Tab. 2 B~ — > 7T~ + 7r° has no asymmetry since the penguins 
cancel due to isospin symmetry. Interesting cases are B~ — > 7r — + 77 and B~ — > 
7T~ + 7/ with branching ratios of the order of 10~ 5 and asymmetries of several percent. 
It is interesting to note that B — ► D~ + D° has a substantOBial rate as well as an 
asymmetry of the order of 2%, independent of whether p is positive or negative. The 
last channel, B" — > X - + fT° is again a pure penguin transition with large asymmetry 
for postive p but with a small branching ratio. 

In Tab. 3 we give the B — > PV results (b — > s transitions). The cases with the 
biggest asymmetries have very small branching ratios (~ 10 -7 ). Interesting decays are 
B~ — > K*~ + 7/ B~ — > K*~ + 77 and B~ — > if* - + 7r° . The three last pairs of decays at 
the bottom of Tab. 3 are again pure penguin modes which have reasonable branching 
ratios but small asymmetries. 

Tab. 4 contains the results for B — > PV (b —>■ d transitions). The last three pairs 
of decays in this table are again pure penguin modes. The transition B~ — ► 7r~ + <p 
has zero asymmetry since no (absorptive) penguin-like matrix elements are present in 
the contributing factorization due to its color structure. The other pure penguin modes 
show appreciable asymmetries for the positive p solution. Their branching ratios are 
small, however. Concerning the other decays, the channel B~ — > p~+r)' stands out with 
a branching ratio of 10 -5 and an asymmetry of 14% for positive p and a branching ratio 
of 6 x 10~ 5 with an asymmetry of 1.2% for the negative p solution. Other interesting 
cases are B~ — > p~ + 77 , B~ — > p~ + ir° , B~ — > n~ + uj and B~ — ► 7r _ + p° with 
significant asymmetries. For p < the asymmetries are usually smaller than for p > 
but then the branching ratios are larger. 

Of course, these results can not be considered definitive since we had to make a series 
of model assumptions namely factorization, BSW matrix elements and the modelling 
of the absorptive contributions in terms of penguin matrix elements. Factorization 
is receiving increasing experimental support |2T|] and also most of the BSW matrix 
elements have been tested experimentally in unsuppressed decays [[2^] . We are aware of 
the colour suppression problem concerning the sign of c\ + C2/N [PT] , [2^| . But we prefer 
to stay in line with the conventional QCD framework together with factorization and 
the limit iV — ► 00. 

As in the case of B — > VV decays, the inclusion of 1/N terms would decrease the 
branching ratios for each of the channels involving a cc- meson by about a factor of 25. 
Other channels which are drastically effected by including 1/N terms are B~ — > 7r~$, 
and B~ — > K~oj, where the branching ratios are decreased by four orders of magnitude 
and a factor of about 20, respectively. In all other channels, the effect of treating the 
1/N terms amounts to a change by less than a factor of two for the prediction of the 
branching ratios or asymmetries (which generally change in the opposite direction). 

A major simplification in our model is that we evaluated the matrix elements at 
fixed k 2 = ml/2. This determines the imaginary parts of the penguin diagrams. In 
a more elaborate treatment the effective k 2 would be determined by the dynamics. 
Here we relied on the work of [26] and |14| yielding a range of values given in fllSD of 
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which we have chosen the upper limit. To get an idea how the asymmetries change 
with k 2 we have calculated the asymmetries as a function of k 2 /m 2 . The results are 
shown in Fig. 2, 3 and 4. Starting at zero for k 2 = the asymmetries vary quite 
appreciably with k 2 (depending on the thresholds appearing in the matrix elements), 
even when k 2 is restricted to the interval (|T%D . Only for the transitions B~ — > D~ + D° 
and B~ — > D~ + D°* the asymmetry is essentially independent of k 2 . It is clear that 
the results presented in the tables in most cases are sensitive to k 2 . A more detailed 
investigation, channel by channel, is needed to get a better estimate of the effective 
value of k 2 , for example, using the approach of Brodsky et al. [13| which includes extra 
gluon exchange or the model of Greub and Wyler [30] using wave functions. 

In Tab. 1, 2, 3 and 4 we have given also results for some Cabibbo allowed channels 
B~ -> K-+r] c , D- + D , D- + D , K- + J/^ } D° + D* s ~, K*~+rj c , D*° + D;, D° + D*- 
and D*° + D~ for completeness. These results can be used for testing some of our 
model assumptions. In these cases the asymmetries are not significant except perhaps 
for B~ — > D° + D*~ and B~ — > D*° + D~ . The channels with two charmed mesons in 
the final state could also be calculated in the frame work of the heavy quark effective 
theory which would be more rigorous than using BSW matrix elements. However, it 
would give similar branching ratios [3jJ and asymmetries (assuming factorization and a 
fixed value of k 2 ). 

We mention that some of the modes have been calculated by other authors using 
a similar approach. First there is work on pure penguin modes |15|, [32], |33| and also 
some work on modes where penguins interfere with tree diagrams @]. Our results are 
in qualitative agreement with these efforts. The work presented here is more elaborate, 
because it includes NLO corrections in the QCD coefficients for the tree diagrams and 
in the penguin contributions which were not taken into account in [34|. Therefore it is 



not too meaningful to make detailed comparisons with this work. The k 2 dependence 
of asymmetries for some selected channels has been studied also |$5], |34|, [15], |32| with 
similar results to ours. 

The branching ratios we have presented in Tab. 1-4 are well within experimental 
limits. However there are three channels where the most recent limits come within 
factors of 2-3 of our predictions: K°tt + , K*°tt + , K + tt°. It is therefore clear that the 
next generation of experiments can provide stringent tests of our model. 
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4 Appendix 

The factorization coefficients z[ 9 l defined in (|T5]) are listed in Tab. [5] for B — > PP 
channels. For the case of B — > PV decays the factorization coefficients are obtained 
by replacing one of the pseudoscalars by a vector meson with the corresponding flavour 
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content, and by setting R[M, P] = if M is a vector meson. 
In Tab. |5] we refer to the unphysical states 

rju = ~^( uU + dd ) , r] s = ss , (20) 

which are related to the physical ^-mesons by taking into account corresponding mixing 
angles. The table covers both cases, q = u and q — c, of ([15]) ; by use of the Kronecker 
delta the contributions from the operators 0\i are to be dropped when tree contribu- 
tions are absent for the given CKM prefactor, see (0). Colour suppressed terms may 
readily be included in the coefficients 

<k = cf + ^cf (21) 

where {i,j} is any of the pairs {1,2}, {3,4}, or {5,6}. In Tab. | we have adopted 
the convention of including factors of V2 associated with a neutral meson Pi. They 
arise either from current matrix elements between P2 and B (left column), or from the 
definition of the decay constants for P2 (right column). Care should be taken with 
the latter since these factors are sometimes absorbed into the decay constants (e.g. as 
tabulated in [0). 
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Figure captions 



Fig. 1: Two types of one- loop matrix elements: (a) Vertex corrections, and (b) penguin 
diagrams. The square box denotes an insertion of one of the four-quark operators 
of §). 

Fig. 2: Momentum dependence of the asymmetry parameter acp for B — > PP transi- 
tions and a p positive CKM Matrix. 

Fig. 3: Momentum dependence of the asymmetry parameter acp for B — > PV (b — > s 
transitions) and a p positive CKM Matrix. 

Fig. 4: Momentum dependence of the asymmetry parameter acp for B —>■ PV (b — > d 
transitions) and a p positive CKM Matrix. 

Table captions 

Tab. 1: Branching ratios and rate asymmetries for B — 
positive and a p negative CKM Matrix, and for k 2 
correspond to the case without strong phases. 

Tab. 2: Branching ratios and rate asymmetries for B — 
positive and a p negative CKM Matrix, and for k 2 
correspond to the case without strong phases. 

Tab. 3: Branching ratios and rate asymmetries for B — 
positive and a p negative CKM Matrix, and for k 2 
correspond to the case without strong phases. 

Tab. 4: Branching ratios and rate asymmetries for B — 
positive and a p negative CKM Matrix, and for k 2 
correspond to the case without strong phases. 

Tab. 5: Factorization coefficients z[ 9 l for various B — > P1P2 decays. The short distance 
coefficients are defined in (|2lD and the factor R is given in (|14|) . 



PP (b — > s transitions) for a p 
= ml/2. Values in parentheses 

PP (b — > d transitions) for a p 
= ml/2. Values in parentheses 

PV (b — > s transitions) for a p 
= ml/2. Values in parentheses 

PV (b — > d transitions) for a p 
= ml/2. Values in parentheses 
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Table 1: 





B — > PP (b - 


s transitions) 




Matrix Elements with (without) Strong 


Phases for k 2 = 


--ml/2 


CKM Matrix: 


p = 0.32,7] 


= 0.31 


p = -0.41, 77 = 0.18 


Channel 


BR 


a C p[%] 


BR 


a C p[%] 


B~ — > K~ + 77 


1.4 xl(T B 


-3.2 


2.0 xlO -7 


-12. 


B+ -> K+ + 7] 


1.5 xl(T 6 




2.5 xl0~ 7 






(1.4 xl(T 6 ) 




(2.2 xl0~ 7 ) 




B~ -> K~ + 7/ 


7.8 xl0~ 6 


5.7 


1.0 xl0~ 5 


2.5 


B+ if + + 7/ 


7.0 xl0~ 6 




9.4 xl0~ 6 






(6.6 xl(T 6 ) 




(8.9 xl0~ 6 ) 




B- K- + 7T° 


6.0 xl0~ 6 


8.1 


1.1 xl0~ 5 


2.3 


B + if + + ^0 


5.1 xl0~ 6 




1.1 xl0~ 5 






(5.3 xl0~ 6 ) 




(1.1 xl0~ 5 ) 




B" K- + Vc 


1.9 xl0~ 3 


0.0 


1.9 xl0~ 3 


0.0 




(1.9 xl0~ 3 ) 




(1.9 xlO- 3 ) 




B~ -> D s + D° 


7.1 xl0~ 3 


-0.12 


7.1 xl0~ 3 


-0.069 


B + -> +W 


7.1 xl0~ 3 




7.1 xlO -3 




(7.1 xl0~ 3 ) 




(7.1 xl0~ 3 ) 




B~ -> ir- +R0 


2.0 xl0~ 5 


0.51 


1.9 xl0~ 5 


0.31 


B + -> 7T+ + K° 


2.0 xl0~ 5 




1.9xl0~ 5 






(1.9 xl0~ 5 ) 




(1.8 xlO^ 5 ) 
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Table 2: 





5 -► PP (6 - 


> d transitions) 




Matrix Elements with (without) Strong Phases for k 2 = 


--ml/2 


CKM Matrix: 


/9 = 0.32,7/ 


= 0.31 


p = -0.41,7? = 0.18 


Channel 


BR 


a C p[%] 


BR 


a C p[%] 


B~ — > 7r~ + 77 


6.6 xl(T B 


-13. 


2.1 xl0~ e 


-21. 


B + -> 7T+ + r? 


8.6 xl0~ 6 




3.2 xl0~ 6 






(7.5 xl(T 6 ) 




(2.4 xl0~ 6 ) 




5" -> + r/' 


1.9 xl0~ 5 


-12. 


3.3 xl0~ 5 


-4.3 


P + -► 7T - + 7/ 


2.5 xl0~ 5 




3.6 xl0~ 5 






(2.2 xl(T 5 ) 




(3.3 xl0~ 5 ) 




5- -► TP - + 7T° 


5.6 xl0~ 6 


0.0 


5.7 xl0~ 6 


0.0 




(5.6 xl(T 6 ) 




(5.7 xl0~ 6 ) 




B~ — »■ 7T~ + 7? c 


8.4 xlO" 5 


0.0 


7.3 xl0~ 6 


0.0 




(8.4 xl0~ 5 ) 




(7.3 xl0~ 6 ) 




B- ^ D- +D° 


4.0 xlO" 4 


2.2 


3.4 xlO~ 4 


1.5 


B+ -> D+ + W 


3.8 xlO~ 4 




3.3 xlO" 4 






(3.9 xlO~ 4 ) 




(3.4xl0- 4 ) 




B- -> K~ + if 


6.8 xl0~ 7 


-17. 


2.3 xlO" 6 


-3.3 


P+ -► + ^ 


9.4 xl0~ 7 




2.5xl0~ 6 






(7.8 xl0~ 7 ) 




(2.2 xlO" 6 ) 
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Table 3: 





B^PV (b- 


s transitions) 




Matrix Elements with (without) Strong Phases for k 2 = 


ml/2 


UK1V1 Matrix: 


p = 0.32, 77 


= 0.31 


p = -0.41, 77 = 0.18 


Channel 


BR 


acp[%\ 


BR 


a C p[%] 


B~ — »■ K~ + u 


7.4 xlO -6 


1.0 


4.2 xlO -6 


1.0 


— ► A ^ + w 


7.3 xlO" 6 




4.1 xlO" 6 






(7.3 xlO" 6 ) 




(4.1 xlO" 6 ) 




73 — > K + p° 


2.7 xlCT 7 


34. 


9.4 xlO" 7 


4.4 


i3 T — > A T + p 


1.4 xlO 




8.6 xlO" 7 






(1.8 xlO" 7 ) 




(8.8 xlO ') 




7~) — T ^ — , Til 

B — ► A + J/^> 


9.0 xlO 


0.0 


9.1 xlO" 4 


0.0 




(9.0 xlO" 4 ) 




(9.1 xlO ) 




5- -► + 79° 


8.8 xlO 


— 0.047 


0.0 X 1U 


—0.02 ' 


E>+ 7~~>*~l~ | 7-1 

IV — ► Dg + l) l> 


8.8 xlO 




8.8 xlO -3 






(8.8 xlO °) 




(8.8 xlO °) 




B — ► iv + 77 


3.4 xlO" 7 


45. 


9.9 x 1U 


6.7 


d — > A +77 


l.o XlU 




8.7 xlO" 7 






(2.2 xlO ) 




(9.2 xlO ') 




7~> — T^H= — 1 / 

B — ► AT + 77 


1.5 xlO" 6 


13. 


2.0 X 1U 


04. 


£> — > A -]- Tj 


1 w 1 n — 6 
1.2 XlU 




5.7 xlO" 8 






(1.3 xlO J 




(9.7 x 10 ) 




7} — > A + 7T 


3.9 XlO 


19. 


1.1 xlO" 5 


3.4 


5+ _^ K*+ + Ti-0 


2.6 xlO" 6 




1.0 xlO" 5 






(3.1 xlO" 6 ) 




(1.1 xlO" 5 




5- -► if*" + 7] c 


5.8 xlO~ 4 


0.0 


5.8 xlO" 4 


0.0 




(5.8 xlO" 4 ) 




(5.8 xlO" 4 ) 




B- -> D s + L>*° 


4.0 xlO" 3 


-0.013 


4.0 xlO" 3 


-0.007 


73+ -► L>J + D*° 


4.0 xlO" 3 




4.0 xlO" 3 




(4.0 xlO" 3 ) 




(4.0 xlO" 3 ) 




B- -> 7T" + if* 


1.5 xlO" 5 


0.54 


1.4 xlO" 5 


0.33 


73 + ->■ 7T+ + K*° 


1.5 xlO" 5 




1.4 xlO" 5 






(1.4 xlO" 5 ) 




(1.3 xlO" 5 ) 




B~ ->K~ + (f) 


5.8 xlO" 6 


1.2 


5.5 xlO" 6 


0.75 


B + -> K+ + (f) 


5.7 xlO" 6 




5.4 xlO" 6 






(4.5 xlO' 6 ) 




(4.3 xlO" 6 ) 




B- -> p- 


1.0 xlO" 6 


0.66 


9.8 xlO" 7 


0.41 


B+ -> p+ + K° 


1.0 xlO" 6 




9.7 xlO" 7 






(9.4 xlO" 7 ) 




(8.9 xlO" 7 ) 
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Table 4: 





B — > PV (b — 


d transitions) 




Matrix Elements with (without) Strong Phases for k = 


m 2 h /2 


CKM Matrix: 


p = 0.32,7? 


= 0.31 


p = -0.41, r 


= 0.18 


Channel 


BR 


a C p[%] 


BR 


a C p[%] 


B~ — > p~ + rj 


7.8 xl(r y 


6.5 


1.8 xl0~ 5 


1.6 


B + -> p+ + 77 


6.8 xl0~ 6 




1.7 xl0~ 5 






(7.3 xl(T 6 ) 




(1.8 xl0~ 5 ) 




5- -► p~ + 77' 


9.6 xl0~ 6 


14. 


5.7 xl0~ 5 


1.2 


B + -► p+ + 7?' 


7.2 xl(T 6 




5.6 xl0~ 5 






(8.2 xl0~ 6 ) 




5.6 xl0~ 5 ) 




B~ ->■ p~ + ^° 


2.3 xlO" 5 


-2.0 


1.7 xl0~ 5 


-1.5 


5+ - p+ + ^° 


2.4 xlO" 5 




1.8 xl0~ 5 






(2.3 xl0~ 5 ) 




(1.8 xl0~ 5 ) 




B~ -► 1T~ + J/Y> 


4.3 xl(T 5 


0.0 


3.7 xl0~ 5 


0.0 




(4.3 xl0~ 5 ) 




(3.7 xl0~ 5 ) 




B~ -► D*- + L>° 


5.0 xlO" 4 


0.89 


4.7 xlO" 4 


0.55 


B+ -> D* + + W 


4.9 xlO" 4 




4.7 xlO" 4 






(5.0 xlO" 4 ) 




(4.7 xlO" 4 ) 




B~ — ► TT~ + W 


1.4 xl0~ 6 


-15. 


1.3 xl0~ 6 


-9.7 


B + -► 7T+ + W 


1.8 xl0~ 6 




1.6 xl0~ 6 






(1.6 xl0~ 6 ) 




(1.4 xl0~ 6 ) 




B~ — > ir~ + p° 


1.3 xl0~ 6 


11. 


2.9 xlO" 6 


2.5 


B + _ 7T+ + p° 


1.0 xl0~ 6 




2.8 xl0~ 6 






(1.2 xl0~ 6 ) 




(2.9 xl0~ 6 ) 




B~ ^ p~ + r) c 


2.7 xl0~ 5 


0.0 


2.4 xl0~ 5 


0.0 




(2.7 xlO" 5 ) 




(2.4 xlO" 5 ) 




Br -> D~ + £>°* 


2.2 xlO" 4 


0.23 


2.2 xlO" 4 


0.14 


B+ -> Z>+ + L>°* 


2.2 xlO" 4 




2.2 xlO" 4 






(2.2 xlO~ 4 




(2.2 xlO~ 4 ) 




B~ — ► 7T~ + 


1.1 xl0~ 7 


0.0 


4.0 xl0~ 7 


0.0 




(1.1 xlO" 7 ) 




(4.0 xlO- 7 ) 




B~ if*- + if 


3.5 xl0~ 8 


-21. 


1.2 xl0~ 7 


-4.3 


B+ -► K* + + ^0 


5.4 xl0~ 8 




1.3 xl0~ 7 






(4.2 xl0~ 8 ) 




(1.1 xl0~ 7 ) 




B- K~ + if * 


4.7 xlO" 7 


-18. 


1.6 xl0~ 6 


-3.5 


B+ -► if+ + if * 


6.7 xlO- 7 




1.7 xl0~ 6 






(5.5 xl0~ 7 ) 




(1.5 xl0~ 6 ) 
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Table 5: 



Pi 


P2 




7(5) 


K~ 
K~ 
K° 
K~ 

Ds 
K~ 


7T° 

Vu 
7T~ 

Vs 

D° 

Vc 


(a 2 S qu + a 4 + a 6 R[K~, tt ])/^ 
(a 2 S gu + a 4 + a 6 i?[K _ ,r? u ])/V2 
a 4 + a 6 R[K°,TT~] 








a\8 qu /\/2 
(aiSqu + 2a 3 + 2a 5 )/\/2 


a 3 + a 4 + a 5 + a 6 J?[t/ s , if - ] 
a 2 5g C + a 4 + a 6 

ai^gc + 03 + 05 


K° 
D~ 

IT" 


7T U 
Vu 

K~ 

Vs 

D° 

Vc 


(fl2<V + a 4 + a6-R[7r _ ,7r°])/\/2 
(a 2 <V + a 4 + a 6 R[tt~,Vu])/ V% 
a 4 + a 6 R[K ,K-} 








(aiS qu - a 4 - a6-R[vr°,7r"])/\/2 
(ai 5qu + 2a 3 + a 4 + 2a 5 + a 6 7r _ ])/V2 


a 3 + a 5 
a 2 <5 9C + a 4 + a 6j R[Z>-, J D ] 
ai5 9C + a 3 + a 5 
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